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1.0  EXECUTIVE  SUMMARY 


This  report  summarizes  the  results  of  a  3-yr  investigation  into  the  role  of  light-scattering  noise  in 
nmliiiear  c^tical  processes.  These  results  span  a  full  range — ^from  physical  and  mathematical 
description  and  prediction  of  noise  properties  to  experimental  verification  of  the  theory  and 
interpretation  of  measurements.  Light-scattering  noise  fluctuations  are  shown  to  set  fundamental 
limits  cm  device  size  and  sensitivity  and  on  beam  power  requirements  for  nonlinear  optical 
processes.  Fundamental  limits  to  optical  phase  conjugation  and  to  weak  signal  amplification  via 
two-wave  mixing  are  derived  through  a  statistical  thermodynamic  treatment  based  on  die 
fluctuation-dissipation  theorem.  In  particular,  power  dissipation  m  a  nonlinear  tptical  medium 
during  the  coherent  transfer  of  energy  from  a  punp  to  a  signal  beam  will  be  examined,  and 
quantitatively  relate  this  dissipation  to  the  sp^trum  of  quantum/thennal  fluctuations  that  give 
rise  to  light-scattering  noise.  Several  classes  of  nonlinear  optical  materials  are  discussed, 
including  resonant  and  nonresonant  Kerr  media  and  photorefractiv'  glasses.  Results  from 
experiments  carried  out  at  the  Science  Center  and  in  collaboration  with  the  University  of  Rome 
(Rmne,  Italy)  ate  analyzed  and  compared  with  computer  simulations  of  beam  propagation  and 
light-scattering  noise  using  the  stochastic  noise  model. 

The  scientific  output  from  the  noise  investigations  includes  two  invited  presentations  (Snowbird 
and  SPIE),  numerous  contributed  papers,  and  many  publications  that  have  appeared  in  journals 
such  as  Physical  Review  (four  papers).  Applied  Physics  Letters  (two  papers)  and  JOSA  B  (two 
piq)ers).  Several  additional  papers  are  in  the  final  stages  of  preparation.  The  project  has  been 
extremely  productive  both  from  our  own  viewpoint  and  from  the  many  independent  reviewers 
who  have  assessed  the  importance  of  our  work,  as  well  as  from  experimentalists  who  have 
undertaken  extensive  measurements  to  confirm  and  extend  our  results.  The  experiments  of 
Stemklar  et  al.  (Ref.  1),  in  particular,  present  a  new  application  of  our  results  to  Brillouin  two- 
beam  coupling  and  continn  our  prediction  of  strong  signal-to-noise  ratio  dependence  on 
pump  intensity.  Eariier  studies  of  Brillouin  amplifiers  (e.g.,  Bespalov  et  al.)  neglected  pump 
dqiendence.  The  experiments  of  Stemklar  et  al.  confirm  our  fluctuation-dissipation  relation  for 
the  amplitude  of  dielectric  fluctuations  in  terms  of  the  temperature,  beam  interaction  volume,  and 
Kerr  coefficient  (Refs.  2-8).  More  detailed  confirmation  can  be  found  in  the  experiments  of 
PizztMfetrato  et  al.  (Ref.  9),  which  also  examine  the  frequency  dependence  of  the  noise 
accon^ianying  two-wave  mixing  in  a  Kerr  medium.  The  experiments  of  Chang  et  al.  (Ref.  1 1) 
extend  the  measurements  of  noise  to  photorefractive  media  and  confirm  the  predictions  obtained 
for  this  class  of  materials  under  the  present  contract. 


1 


It  is  believed  that  with  further  development  the  results  will  serve  as  the  scientific  basis  necessary 
for  the  (tesign  of  new  low-noise  nonlinear  optical  materials — ^probably  with  best  prospects  in  the 
general  classes  of  photorefractive  polymers  and  charge  transfer  media.  Furthermore,  our  results 
can  provide  the  foundation  for  future  investigations  into  the  noise  properties  of  optical 
communication  channels  and  optical  storage  media.  These  future  investigations  can  include 
determination  of  the  Shannon  channel  capacity  (some  work  along  these  lines  has  already  been 
carried  out  for  photorefractive  media)  and  new  applications  of  optical  signal  processing  and 
ctxnputing  (e.g.,  the  optical  wavelet  transform  and  optical  methods  for  signal  and  image- 
conqnession). 
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2.0  INTRODUCTION 


Most  treatments  of  (q>tical  phase  conjugation  and  two-wave  mixing  do  not  address  the  role  of 
noise;  others  include  quantum  noise  in  the  incoming  electromagnetic  fields,  but  neglect 
temperature-dependent  noise  arising  fo»n  fluctuations  inherent  in  the  nonlinear  medium  itself. 
Such  fliKtuaticnis  are  itx^xnrtant  to  consider  to  the  extent  that  they  reduce  the  fidelity  of  an  optical 
signal  by  giving  rise  to  scattered  light. 

• 

Studies  were  done  on  light-scattering  noise  for  continuous  wave  (cw)  nonlinear  optical  processes 
using  Kerr  and  artificial  Kerr  media  (Refs.  2-9)  and  for  similar  processes  using  the 
photorefractive  effect  (Refs.  10  and  1 1).  In  the  former  case,  a  fundamental  connection  between 
light  scattering  fluctuations  Se  in  the  medium  and  the  nonlinear  dielectric  constant  £2  was 
derived  in  the  form  of  a  static  susceptibility  relation  (Ref.  4)  and  was  applied  to  a  study  of  light¬ 
scattering  noise  in  the  phase  conjugate  signal  obtained  through  degenerate  four-wave  mixing 
(Refs.  5  and  6).  Time-averaged  noise  powers  were  determined  and  were  of  order  kTv  where  k  is 
the  Boltzmann  constant,  T  is  temperature,  and  v  is  the  optical  fiequency.  At  room  temperature 
and  visible  wavelengths,  kTv  is  in  the  microwatt  range.  The  static  susceptibility  between  £2  and 
5e  was  extended  in  a  more  recent  article  and  applied  to  four-wave  mixing  in  Kerr  media  with 
nonlocal  interactions  near  a  critical  point  (Ref.  S).  Throughout  the  following,  the  noise  arising 
from  thejfluctuations  5£  equivalently  will  be  referred  to  as  either  thermal  or  light-scattering 
noise.  • 


For  many  optical  signal  processing  applications  a  more  complete  description  of  noise  than  is 
furnished  by  the  time-averaged  noise  power  is  desired.  This  additional  information  might 
include,  for  example,  the  effect  of  light-scattering  noise  fluctuations  on  the  output  light  intensity 
of  a  nonlinear  c^tical  device  in  real  time.  For  image  processing  applications,  such  as  the  use  of 
optical  phase  conjugation  to  correct  aberration,  the  effects  of  light-scattering  noise  on  the 
anq)litude  and  phase  of  the  electromagnetic  field  itself  must  be  included.  In  the  following 
sections,  stochastic  simulations  for  light-scattering  noise  are  developed;  the  simulations  provide 
a  direct  means  through  which  this  addititmal  infrmnation  may  be  obtained.  In  the  following 
sections,  the  theory  of  Ught-scatteiing  noise  in  Kerr  media  is  extended  using  the  fluctuation- 
dissipation  theorem  (Ref.  6).  This  generalization  of  the  static  susceptibility  relation  is  required 
for  triplications  to  time-dependent  processes  and  nonlinear  optical  media  exhibiting  power 
dissipation.  An  important  example  of  such  a  process  is  the  amplification  of  a  weak  optical  signal 
via  die  coherent  transfer  of  energy  from  a  strong  pump  beam  during  nondegenorate  two- wave 
mixing  (Refs.  7-9).  Reference  6  presents  the  basis  for  stochastic  noise  simulation  using  Langevin 
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and  generalized  Langevin  models  to  relate  the  time  decay  of  laser-induced  gratings  to  the  thermal 
flucmations  lesptmsible  for  light-scattering  noise.  Applications  to  the  nonlinear  optical  processes 
of  phase  conjugation,  via  degenerate  four-wave  mixing  (Ref.  6),  and  to  weak  signal 
anq)lification,  via  nondegenerate  two-wave  mixing  (Refs.  7-9),  were  obtained  as  part  of  the 
study. 

Figure  1  shows  a  schematic  depiction  of  the  origins  of  light-scattering  noise.  The  solid  arrows 
show  incoming  and  outgoing  signal  beams  with  a  volume  element  of  the  nonlinear  medium 
while  the  dashed  arrows  show  the  scattered  light  noise.  The  amplitudes  and  time-correlation 
properties  of  the  signal  and  noise  gratings  (fluctuations  in  the  dielectric  constant  of  the  medium) 
are  connected  through  the  fluctuation-dissipation  relation  described  below.  For  Kerr  media,  the 
signal  response  and  noise  properties  are  connected  via  the  fluctuation-dissipation  theorem  to 
thermal  fluctuations  in  density.  For  photorefractive  media,  the  signal  is  connected  to  thermal 
fluctuations  in  the  space-charge  field.  Photorefractive  noise,  however,  is  due  to  scattered  light 
from  both  thermal  fluctuations  in  the  space-charge  field  and  thermal  fluctuations  in  density.  Thus 
for  a  photorefractive  medium,  the  Kerr  coefficient  needs  to  be  determined  to  completely  describe 
the  noise  (Refs.  10  and  1 1). 


Photortfrictlve  -  Thtrmal  Fluctuations 

In  tlia  Spaca-Charge 
Field 


Karr  Effect  - 

% 

— - ■*  Thermal  Fluctuations 

In  Density 

Figure  1.  Light-scattering  noise. 
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3.0  LIGHT-SCATTERING  NOISE  IN  PHOTOREFRACTIVE  MEDIA 


Photorefractive  media  constitute  a  large  and  versatile  class  of  nonlinear  optical  materials  with 
applications  to  holographic  storage  as  well  as  to  real-tin»  optical  information  processing.  Two 
closely  related  theoretical  models  of  the  photorefractive  effect  have  been  put  forth  (Refs.  12  and 
13)  and  successfully  used  over  the  past  decade  by  a  number  of  investigators  for  the  description  of 
nonlinear  optical  processes  including  optical  phase  conjugation  and  coherent  beam  combination 
via  two-wave  mixing.  Nevertheless,  few  studies  have  been  directed  at  quantifying  the 
fundamental  limits  to  the  photorefractive  effect  set  by  noise.  In  particular,  questions  related  to  the 
noise  limits  on  photorefractive  sensitivity  to  weak  incident  signal  powers  (dynamic  range) 
apparently  have  not  been  addressed.  The  present  study  shows  that  noise  limits  to  the  dynamic 
range  of  a  photorefractive  medium  are  inherent  in  the  photorefractive  effect  itself.  The 
quantitative  determinadon  of  these  limits  is  the  foremost  objective  of  this  secdon. 

3.1  ORIGINS  OF  THE  PHOTOREFRACTIVE  EFFECT 

Inidadon  of  the  photorefracdve  effect  occurs  when  a  spatial  variation  of  light  intensity  incident 
on  the  medium  causes  a  redistribution  of  charge  density  and  buildup  of  a  space-charge  field.  For 
example,  let  the  spatial  light  intensity  distribution  follow  the  simple  sinusoidal  grating  form 

I(x)  =  Io(l +mcos(qx)]  (1) 

characterized  by  average  intensity  Iq,  grating  wave  vector  q,  and  intensity  modulation  ratio  m. 
The  intensity  grating  induces  a  corresponding  (tensity  of  the  form 

n(x)  =  noli  +  w  cos(qx)]  (2) 

where  no  is  the  average  (tensity.  In  the  hopping  model,  n  refers  to  the  number  density  of  hopping 
(^Triers  of  charge  q^.  In  the  Kukhtarev  model,  n  refers  to  the  number  density  of  unionized  donor 
sites. 

A  comparison  of  the  Kukhtarev  and  hopping  mcxtels  is  presented  in  Figure  2.  Each  nKxtel 
assumes  a  uniform  population  of  fixed  sites  as  well  as  a  uniform  distribution  of  fixed 
countercharges  to  maintain  overall  charge  neutrality.  In  the  hopping  model,  the  fixed  sites  may 
be  occupied  by  a  smaller  number  of  mobile  carriers  of  charge  qe.  The  latter  are  assumed  to  hop 
firom  site  to  site  with  a  hopping  rate  that  is  proportional  to  the  local  light  intensity.  The  net  effect 
is  to  cause  the  carriers  to  diffuse  from  regions  of  higher  to  regions  of  lower  light  intensity.  As  a 
result,  the  density  modulation  ratio  w  in  Equation  2  is  opposite  in  sign  to  the  intensity 
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modulation  ratio,  m.  In  the  Kukhtaitv  model,  the  fixed  sites  are  occupied  with  a  smaller  number 
of  valence  band  carriers  whose  charge  may  be  either  positive  or  negative.  For  definiteness,  the 
carriers  in  both  models  are  assumed  to  be  electrons  of  charge  =  e  =  -  lei.  Mobility  of  carriers  in 
the  Kukhtarev  model  requires  excitation  from  the  valence  band  into  the  conduction  band 
(depicted  by  the  shaded  region  in  Fig.  2)  at  a  rate  that  is  proportional  to  the  light  intensity. 
Generation  and  recombination  of  carriers  is  a  fast  process  compared  to  the  overall  redistribution 
of  charge,  and  the  number  of  carriers  in  the  conduction  band  at  any  given  time  is  negligible 
compared  to  the  number  of  valence  band  carriers  (donors).  Redistribution  of  charge  in  the 
Kukhtarev  model  may,  therefore,  be  thought  of  as  due  to  the  redistribution  of  the  donor 
population  through  accumulated  generation,  recombination,  and  transport  of  carriers  via  the 
conduction  band.  The  index  change  induced  by  the  space-charge  Held  occurs  through  the  Pockels 
effect  as  described  below. 


o  o 

TTTfT  fTfTT  TTTTT 


-CARRIER  (q) 


rmr  rnn - FIXED  SITES 


FIXED  COUNTERCHARQE 


HOPPING  MODEL 
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COM>UCTION 

BAND 


CARRIER  (•) 
FIXED  SITES 


FIXED  COUNTERCHAROE 


KUKHTAREV  MODEL 


Figure  2.  Photorefractive  glasses:  comparison  of  hopping  and  Kukhtarev  models. 

To  simplify  the  present  analysis,  contributions  to  grating  formation  from  the  photocunent,  dark 
current,  and  externally  applied  electric  rields  will  be  neglected.  The  contributions  from  higher 
order  gratings  will  also  be  neglected  since  the  primary  interest  is  in  the  small  modulation  ratio 
limit  The  various  components  of  the  grating  free-energy  will  now  be  derived.  These  include 
contributions  from  energy  storage  in  the  space-charge  Held,  configurational  entropy,  and  the  free 
energy  associated  with  exposure  of  the  donor  or  hopping  sites,  depending  on  which  model  is 
used,  to  a  spatially  varying  light  intensity. 
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Fluctuations  in  the  space-charge  field  are  determined  by  the  free-energy  change  that  accompanies 
the  fluctuation.  The  free  energy  is,  in  turn,  made  up  of  energy  and  entropy  contributions  as 
described  in  Reference  10.  In  the  diffusion  limit,  the  energy  contribution  dominates  and  thermal 
fluctuations  in  the  space-charge  field  may  be  obtained  immediately  from  the  equipartition  of 
energy  stored  in  the  space-charge  field 

Use  =  (eo/8Jt)  /  Esc^  (x)  dx  =  <  I  5Escl2>  =  kT/2  (3) 

In  the  absence  of  an  entropic  contribution,  Ujc  is  equal  to  the  free  energy.  Rewriting  the  last 
equality  gives 

<  1 5E5g(q)l^>^f  limit  =  (47ikTyiEQV  s)  (4) 

The  general  result  is  (Ref.  10) 

<  1 5Esc(q)|2>  =  (4jckT/feoVs)  [1  +  (qAcD)^]-!  (5) 

where  q  is  the  magnitude  of  the  grating  wave  vector  and  kp  is  the  reciprocal  Debye  length.  The 
difference  between  Equations  4  and  5  is  due  to  the  entropic  contribution  with  the  square  bracket 
factor  accounting  for  nonlocal  interaction  (Ref.  10). 

3.3  LIGHT-SCATTERING  FLUCTUATIONS  ASSOCIATED  WITH  THE 

PHQTQREFRACTIVE  EFFECT 

Dielectric  fluctuations  that  originate  from  fluctuations  in  the  space-charge  field  follow 
immediately  from  Equation  5  through  the  Pockels  effect 

5e  =  e(2)5Esc  (6) 

where  e(2)  is  related  to  the  electro-optic  tensor,  r,  and  the  refractive  index  of  the  medium,  ur.  The 
standard  relations  An  =  -(l/2)nR^E,  for  the  index  change,  and  Ae  =  2nRAn  yield 

e(2)=:-nR4r  (7) 

(jombining  Equation  5  for  the  mean-square  fluctuations  in  the  space-charge  freld  with  Equation  6 
gives 

<  l5e(q)|2>pR  =  (4jckT/eoVs)  [1  +  (q/kD)2]-l[e(2)]2  (8) 
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A  subscript  PR  has  been  added  to  designate  those  fluctuations  associated  with  the  photorefractive 
effect  Equation  8  is  the  central  result  of  Reference  10.  It  provides  the  extension  of  our  previous 
investigations  of  light-scattering  noise  in  Kerr  media  to  materials  exhibiting  the  photorefractive 
effect 

3.4  LIGHT-SCATTERING  FLUCTUATIONS  IN  PHOTOREFRACTIVE  MEDIA 
ASSOCIATED  WITH  THE  OPTICAL  KERR  EFFECT 

Dielectric  fluctuations  arise  in  all  media  (including  photorefractives)  that  exhibit  a  Kerr^ffect 
(Ref.  4)  and,  therefore,  contribute  to  light-scattering  noise.  These  fluctuations  originate  from 
thermal  fluctuations  in  the  density,  which  through  the  Clausius-Mossotti  constant  de/dp  induce 
corresponding  fluctuations  in  the  dielectric  constant  independent  of  the  photorefractive  effect 
(Ref.  5).  From  Reference  5 

<  1 8e(q)|2>K  =  87ckTe2(q)A^s  =  (8jtkTe2A^s)  [1  +  (q/qo)^!-!  (9) 

where  e2(q)  is  the  Kerr  coefficient  and  qo  is  the  reciprocal  correlation  length.  Generally,  for  Kerr 
media  q/qo  «  1  and  this  is  true  for  the  photorefractive  media  whose  noise  properties  were 
studied  (Ref.  11).  The  inequality  breaks  down  when  the  correlation  length  approaches  the 
wavelength  of  light  as  can  occur  close  to  a  critical  point  (the  case  for  which  Eq.  9  was  originally 
derived  in  Ref.  5).  The  purpose  of  including  the  general  result  here  is  to  note  its  similarity  with 
Equation  8  for  the  dielectric  fluctuations  associated  with  the  photorefractive  effect.  Apart  from 
the  physical  distinction  between  the  correlation  length  and  the  Debye  screening  length  in  the  two 
models,  the  forms  differ  only  by  a  nondimensional  factor  2eq  e2/[e^2)]2 

3.5  ESriMATESGEDYNAMIC  RANGE  AND  COMPARISON  WITH  EXPERIMENT 

The  methods  developed  in  the  preceding  subsections  are  illustrated  in  this  subsection  for  the 
write  and  readout  of  a  simple  sinusoidal  grating  using  the  configuration  shown  in  Figure  3.  For 
the  writing  beam  configuration  of  Figure  3,  the  intensity  modulation  ratio  is: 

m  =  2ail2)*'2/(ij  +  i2)  (10) 

where  Ii  and  I2  are  the  intensities  for  the  pump  (Ii)  and  signal  (I2)  beams.  For  light  of  wave 
vector  Ki  =  K2  =  K,  the  resulting  grating  is  of  wave  vector  q  =  Iql: 


q  =  2Ksin(e/2) 


(11) 


where  6  is  the  angle  between  the  writing  beams.  Both  K  and  6  are  measured  in  the  material  with 
refractive  index  nj^. 


Writing  Beam  Ratio  (P) 

Figure  3.  Measurements  of  dynamic  range  (from  Ref.  1 1). 


In  the  read  stage  (Fig.  3),  the  intensity  of  the  diffracted  beam  (signal  plus  noise)  is  proportional 


to  the  square  of  the  total  grating  amplitude 

Ijiff- (aj  +  6ai)2  (12) 

where  aj  and  5a  i  are  the  amplitudes  of  the  signal  and  noise  grating  components,  respectively 
(Ref.  10): 

ai  =  -(kTq/e)  ( 1  +  (q/kD)2]-lE(2)m  (13) 

The  noise  grating  component  is  Gaussian  distributed  with  mean 

<5a|>  =  0  (14a) 

and  variance 

<  1 5ail2>  =  2  <  1 6el2>  (14b) 
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The  factor  of  2  in  Equation  14b  corrects  for  averaging  the  square  of  the  sinusoidal  spatial 
dependence  of  5e  over  the  beam  interaction  volume.  The  right-hand  side  of  Equation  14b  is 
evaluated  using  Equation  8  or  9,  depending  on  whether  the  noise  in  question  is  associated  with 
the  Kerr  or  photorefractive  effects. 

Experimentally,  it  is  expedient  to  conduct  measurements  under  the  condition  that  q  =  kp  since  it 
follows  from  Equation  13  that  the  diffraction  efficiency  is  maximized  when  this  condition  is 
satisried.  Under  this  condition  the  signal-to-noise  ratio  (Ref.  10)  from  is  ngYsm^M,  which  equals 
unity  for  a  writing  beam  intensity  ratio  I2/I1  l/(noVs).  Figure  4  from  Reference  10  shows  a 
calculation  of  light-scattering  noise  in  the  difOacted  read  beam  for  BaTi03  under  the  maximum 
diffraction  efficiency  condition.  The  figure  shows  the  logarithm  of  the  diffiacted  read  beam 
powo*  (in  arbitrary  units  due  to  the  unknown  prefactor  in  Eq.  12)  versus  the  logarithm  of  the 
writing  beam  ratio  I2/I1.  The  beam  interaction  volume  was  set  equal  to  the  product  of  the  beam 
cross  section  area  and  the  interaction  length.  For  the  present  calculation,  the  beam  diameter  was 
set  at  O.IS  cm  and  the  interaction  length  at  0.2  cm,  corresponding  to  the  experimental 
measurement  conditions  used  in  Reference  11.  Other  parameters  used  in  the  calculation,  which 
were  chosen  as  representative  of  BaTi03.  are  listed  in  Reference  10.  As  described  below,  the 
light-scattering  fluctuations  associated  with  the  Kerr  effect  are  relatively  less  important  for 
BaT103  and  were  omitted  from  the  calculation  of  Figure  4.  This  will  not  always  be  the  case  for 
other  photorefractive  materials  as  shown  below. 
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To  obtain  Hgure  4,  the  right-hand  side  of  Equation  14b  for  the  noise  variance  was  evaluated 
using  Equation  8  with  q  s  ki).  To  generate  the  figure,  the  x-coordinate  Gog  writing  beam  ratio) 
was  randranly  selected,  with  unifcmn  distribution  along  the  scale  of  the  figure.  Next,  the  signal 
grating  an^litude  a^  was  obtained  from  Equation  13  with  m  from  Equation  10.  Finally,  a  value 
for  Saj  was  sampled  from  the  Gaussian  distribution  having  mean  and  variance  specified  by 
Equatimis  14a  and  14b.  Gaussian  probability  sampling  was  achieved  using  a  standard  computer 
program  incorporating  the  Box-Muller  transformation  for  the  generation  of  normal  deviates  from 
random  numbers  sampled  uniformly  on  the  interval  (0,1 }.  The  above  procedure  was  repeated  a 
thousand  times  to  obtain  an  equal  number  of  points  for  good  statistical  sampling  of  the  noise. 

Expansion  of  the  right-hand  side  of  Equation  12  gives  three  terms,  each  corresponding  to  a 
different  region  seen  in  Figure  4.  The  first  term  (aj^)  dominates  in  the  high  probe  power  Gow 

ncnse)  region.  In  this  case,  the  linear  behavior  (with  unit  slope)  predicted  from  Equation  10  for 
I2  <  0.1 1)  is  observed.  As  the  probe  power  is  reduced  further,  the  next-to-leading  term  in  the 

expansion  (2ai5ai)  begins  to  play  a  role.  This  term,  which  vanishes  on  ti-'*  averaging,  results  in 
a  symmetric  distribution  of  scatter  about  the  average  diffracted  beam  power  due  to  noise.  This 
behavior  is  seen  in  the  mild  noise  region  of  Figure  4.  With  still  further  reduction  of  probe  power 
(i.e.,  reduction  of  ai).  the  third  term  in  the  expansion  (Saj^)  begins  to  dominate.  This  last  term, 
which  does  not  vanish  on  time  averaging,  results  in  a  diffracted  power  due  entirely  to  noise.  The 
resulting  scatter  is  independent  of  the  power  in  the  probe  since  in  this  region  the  signal  grating 
amplitude  is  much  smaller  than  the  RMS  fluctuation  amplitude  due  to  noise.  This  effect  is  best 
seen  in  the  smallest  writing  beam  ratio  range  of  Figure  4.  The  probe  power  level  for  which  the 
signal-to-noise  ratio  is  unity  may  be  computed  analytically.  For  the  present  values  of  Vs  and  oq, 
(Mie  obtains  I2/I1  =  5.9  x  lO-iS  for  S/NpR  =  1. 

Equations  8  and  9  are  compared  for  their  relative  importance  in  limiting  the  dynamic  range  of  a 
photorefractive  medium.  Setting  (mi^)^  =  a(mj2)pR  one  obtains  for  a: 

a  =  {2£oe2/[e(2)]2}(i+(qAcD)2) 

=  {4e()n2/(nR7r2)}  [  1  +  iq/k^H^]  (15) 

For  BaTiC)3,  the  parameters  collected  in  Reference  10  yield  a  value  of  0.035  for  the  expression  in 
curly  brackets  and  dynamic  range  is  limited  by  thermal  fluctuations  in  the  space-chaige  field.  For 
different  materials  and/or  different  writing  beam  configurations,  the  curly  bracketed  expression 
in  Equation  15  may  easily  exceed  unity;  in  which  case  light-scattering  noise  associated  with  the 
Keir  effect  will  dominate.  Values  of  a  are  computed  for  several  other  photorefractive  materials 
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in  Reference  11.  There  one  finds  a  values  of  0.93  for  KNb03, 1.56  for  BSKNN,  and  3.75  for 
SBN;  all  for  the  maximum  diffraction  efficiency  condition  q/k£)=l  (Table  1).  Thus  for  KNb03 
and  BSKNN,  it  is  predicted  that  the  thermal  light-scattering  fluctuations  associated  with  the  Kerr 
and  photorefiractive  effects  are  of  comparable  size,  while  for  SBN,  fluctuations  associated  with 
the  Kerr  effect  are  of  greater  importance  to  limiting  dynamic  range. 

Table  1.  C^culated  fundamental  limit  on  the  dynamic  range 
due  to  the  various  noise  sources  discussed. 


Type  of 
Noise 

Minimum 
Beam  Ratio 

BaTi03 

BSKNN- 11 

KNb03 

SBN;60 

Photon  Shot 
Noise 

P^.inin 

5.3  X  10-18 

3.1  X  10-17 

4.9  X  10-17 

3.4  X  10-17 

Photorefractive  a _ 

Noise 

2.6  X  10-15 

2.3  X  10-14 

1.4  X  10-15 

6.4  X  10-16 

Kerr  Noise 

Pketr,  min 

3.0  X  10-16 

3.6x10-14 

1.3  X  10-15 

2.4  X  10-15 

Finally,  note  that  there  will  be  a  photon  shot  noise  contribution  arising  from  fluctuations  in  the 
incident  intensity  of  the  weak  writing  beam  when  the  latter  becomes  sufficiently  small  that  few 
photons  from  this  beam  arrive  during  the  time  required  for  photorefiractive  grating  formation 
(Ref.  1 1).  The  photon  shot  noise  limits  to  the  minimum  writing  beam  ratio  were  evaluated  for 
each  of  the  various  photorefractive  materials  listed  above  and  were  in  the  10*1^  to  10*1'^  range 
(Ref.  11)  Thus,  for  these  materials,  the  photon  shot  noise  is  predicted  to  be  about  1-2  orders  of 
magnitude  smaller  than  the  thermal  noise  associated  with  either  the  Kerr  or  photorefractive 
effects.  The  results  presented  in  this  section  suggest  a  very  large  dynamical  range  for 
phtnorefiactive  materials  that  should  prove  useful  in  optical  signal  processing  applications. 

Experimental  measurements  for  the  dynamic  range  using  BaTi03  as  the  recording  medium  are 
shown  in  Figure  3  from  Reference  1 1  (note  the  similarity  to  the  theoretical  calculations  of 
Fig.  4).  In  tills  case  the  measured  dynamic  range  is  large  but  not  as  large  as  predicted  due  to  the 
limitations  of  detector  noise  in  the  experiments.  The  presently  available  measurements,  thereftne, 
serve  mdy  as  an  upper  bound  to  the  fundamental  noise  limits  of  photorefractive  media.  Note  that 
the  noise  properties  of  photorefractives  are  qualitatively  different  from  those  of  Kerr  media.  For 
the  latto’,  light-scattering  noise  is  significant  at  room  temperatures  even  for  signal  powers  in  the 
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milliwatt  range.  Differences  are  predominantly  in  the  physical  mechanisms  underlying  formation 
of  the  signal  grating  (Ref.  10).  Finally,  note  that  analytic  expressions  for  the  signal-to-noise  ratio, 
based  mi  the  preceding  results,  provide  the  means  for  predicting  the  channel  capacity  of  a 
photorefractive  medium  (bits/sec)  for  signal  processing  applications. 
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4.0  GENERAL  THEORY  OF  LIGHT-SCATTERING  NOISE  IN  KERR  MEDIA 


An  overview  of  light-scattering  fluctuations  and  their  connection  to  nonlinear  optical  response 
may  be  gained  through  an  examination  of  the  Maxwell  equations  governing  beam  propagation  in 
a  mmlinear  medium: 


V  X  E(r,t)  =  -( l/c)9H(r,t)/3t  ( 1 6a) 

V  X  H(r,t)  s=  (l/c)3D(r,t)/3t  (16b) 


with 

D(r,t)  =  [  Co  +  8e(r,t)  +  £2!^  (r,t)]  E(r,t)  (16c) 

In  these  equations.  E(H)  is  the  total  electric  (magnetic)  field  vector;  eq  (£2)  is  the  linear 
(nonlinear)  dielectric  constant  of  the  medium;  D(r,t)  is  the  displacement  vector;  and  the  overbar 
implies  a  time  average  that  is  long  compared  to  an  optical  period,  but  shon  compared  to  all  other 
time  scales  that  enter  the  problem.  Equation  16c  contains  both  the  nonlinear  contribution  and  the 
fluctuation  contribution  to  the  dielectric  constant  The  symbols  &£  and  A£  »  are  used 
throughout  this  repart  to  designate  variations  in  the  dielectric  constant  of  the  medium  due  to 
spontaneous  fluctuadons  and  to  field-induced  response,  respectively.  The  fluctuation 
contribution  plays  a  role  in  nonlinear  optics  analogous  to  the  current  or  voltage  fluctuations  that 
give  rise  to  Johnson  noise  in  electrical  systems.  Both  sources  of  noise  are  thermal  in  origin  and 
prcqxntional  to  kT.  In  the  present  cox,  fluctuations  in  the  linear  dielectric  constant  give  rise  to 
noise  as  they  result  in  the  formation  of  spontaneous  gratings  that  are  “read”  by  the  applied  laser 
fields  to  give  a  scattered  light  component  that  is  indistinguishable  from  the  desired  output  of  the 
nonlinear  t^tical  device. 

4.1  LIGHT-SCATTERING  FLUCTUATIONS  AND  NONLINEAR  OPTICAL 

RESPONSE 

At  nonzero  temperature  T.  thermal  fluctuations  in  the  linear  dielectric  constant  of  the  optical 
Kerr  medium  give  rise  to  fluctuation  gratings  capable  of  scattering  incident  radiation.  Mtne 
precisely,  the  fluctuations  in  the  linear  dielectric  constant,  6£(r,t),  can  be  decomposed  into 
grating  components,  5£(q,t): 

5e(r,t)  =  Z  5£(q,t)  e-ifl***  (17) 


with 
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8e(q)  -  ( Ws)J  ciqT8E(r)d3r  (18) 

where  the  summation  is  over  all  grating  wave  vectors  q,  and  the  integral  is  over  the  beam 
interaction  volume  Vs. 

For  an  is(Krq)ic  Kerr  medium  with  local  interactions 

e(r)  =  eo  +  e2E^(r)  ^  (19) 

and  Fourier  transformation  (Eq.  18)  gives 

e(q)  *  eo5(q)  +  £2  E2(q)  (20) 

A  nonlocal  generalization  of  Equation  19  may  be  written  in  terms  of  the  convolution  integral: 

e(r) »  eo  +  /  e2(r  -  r’)E2(r*)  d3r’  (21) 

far  an  isotn^ic  medium.  Fourier  transformation  gives 

e(q)  *  €b8(q)  +  e2iq)EHq)  (22) 

This  result  is  identical  in  fcmn  to  Equation  20  ex<xpt  for  the  q*dependence  in  £2,  which  becomes 
inqKxtant  when  the  range  of  correlations  in  the  medium  approaches  the  wavelength  of  light 
(Ref.  5). 

The  remainder  of  this  section  shows  that  the  fluctuation-dissipation  theorem  provides  a 
fundamental  closure  relation  between  the  nonlinear  coefficient  ex  and  the  statistical  distribution 
determining  the  fluctuations  Se  for  incorporation  into  Equation  16c. 


Preliminary  to  ^^lying  the  fluctuation-dissipation  theorem,  the  appropriate  conjugate  variables 
must  be  defined.  For  this  purpose,  it  is  sufficient  to  note  that  the  energy  change  in  the  medium 
due  to  the  nonlinear  polarization  is  of  the  form 

U(r.t)  *  .(l/8JC)Ae(r,t)  E^  (r,t)  (23) 

whoe  Ae  is  the  change  in  the  dielectric  constant  in  the  presence  of  an  applied  laser  field. 
Accodingly,  the  field  energy  variable  u(r,t)  -  (l/8ic)  E^  is  conjugate  to  Ae  and  one  may  define 
the  conqilex  frequency-dependent  susceptibility  %£  such  that: 
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Ac(q,Q)  =Xe(Q.ii) 


(24) 


witfi  Xe  *  Xe'  ^  i  Xe"-  particular,  it  follows  ftom  Equation  24  for  a  Kerr  medium 
(Ae(q,0) »  62(q)E2(q))  that  XgCq.O)  =  8jc  ejfq)  =  Xe'(q*0)- 

The  power,  P,  expantted  by  the  applied  fields  on  the  medium  electric  polarization  is  proportional 
to  the  imaginary  component  of  the  susceptibility 

p  =  <E*dPNL/dt>  =  VsftXe"u^/2  *  (25) 

where  Pjg^  ts  the  nonlinear  polarization  and  the  angular  brackets  denote  averaging  over  the  beam 
interaction  volume  V^.  The  fluctuation-dissipation  theorem  is  (Ref.  6) 

Xe"  (q.Q)  =  (l/2h)  [  1  -  exp(-hft/kT) )  S^iqjO)  (26a) 

which  in  the  classical  limit  (kT  »  hQ),  reduces  to  the  form 

Xe"(q^)  =  (a/2kT)  Se(q,Q)  (26b) 

v^ieie 

Se(q.O)  =  Vs  J  <8e(q,0)  8e(q,t)  >  dt  (27) 

is  the  spectral  density  of  the  fluctuations  in  the  linear  dielectric  constant  5e  that  give  rise  to  light 
scattering  noise. 

Each  of  die  preceding  results  may  readily  be  expressed  in  tensor  component  form  for 
aiqilkations  to  anisotropic  media.  The  scalar  u,  for  example,  is  then  replaced  by  the  Maxwell 
stress  tensor  in  its  most  general  form,  with  ^  and  Xe  represented  by  second  and  fourth  ranked 
tensors,  respectively. 

Equations  25  through  27  are  the  appropriate  generalization  of  our  previously  derived  static 
susceptibility  relation,  which  was  limited  to  degenerate  beam  interactions  and  did  not  include 
powOT  dissipation.  To  recover  the  special  case,  first  note  that  the  real  and  imaginary  components 
of  Xe  satisfy  the  Kramers-Kronig  relation  (Ref.  6): 

Xe‘  (q.O)  =  (IM)/  Xe"  (q^)  /«  ^ 

=  (kT)-lVs<l6e(q)|2> 
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(28) 


where  the  last  equality  follows  upon  substitution  from  Equations  26b  and  2^  for  the  imaginary 
susceptibility  component  and  the  integral  is  over  both  positive  and  negative  frequencies.  Rom 
Equation  28,  and  the  above  assignment  of  Xe(4>0)  for  a  Kerr  medium, 

Xe'(<I*0)*8jce2(q)  (29) 

the  static  susceptibility  relation  (Ref.  4) 

<  I8e(q)|2>  =  87ckTe2(q)/Vs  .  (30) 

is  immediately  obtained.  Equation  30  is  the  closure  relation  required  for  incorporation  into 
Equatitm  16c.  Its  present  derivation  is  both  more  direct  and  more  general  than  the 
thermodynamic  fluctuation  approach  used  in  Reference  S,  which  also  allows  for  possible 
q-dependence  in  the  Kerr  coefficient  due  to  nonlocal  interaction. 

For  the  case  of  photorefractive  media  described  in  the  preceding  subsection,  dielectric 
fluctuatitms  arise  from  thermal  fluctuations  in  the  space-charge  field  through  the  electro-optic 
effect  (Ref.  10).  In  this  case,  fluctuation  variance  is  given  by  an  expression  which  differs  from 
the  right-hand  side  of  Equation  30  only  by  the  nondimensional  factor  [e(2)]2/2eQe2(0)  where  e(2) 
is  the  Pockels  coefficient  and  the  q  dependence  takes  the  same  form  as  that  obtained  previously 
for  the  nonlocal  interaction  described  in  Reference  5,  with  the  range  of  correlations  set  by  the 
Debye  smeening  length. 


5.0  APPLICATIONS  TO  TWO-WAVE  MIXING  AND  COMPARISON  WITH 

EXPERIMENT 


The  first  measurements  of  noise  during  two-wave  mixing  in  a  Kerr  medium  are  reported  in 
References  8  and  9.  Experiments  were  conducted  at  visible  wavelengths  using  an  aqueous 
suspension  of  shaped  microparticles  as  the  nonlinear  medium.  Theoretical  calculations  and 
C(xiq)uter  simulations  using  a  stochastic  model  for  light-scattering  noise  gave  excellent 
agreement  with  the  experimental  results.  This  section  presents  a  summary  of  those  resufts. 

5.1  IHEQRY 

A  schematic  illustration  of  the  two-wave  mixing  geometry  for  amplification  of  a  weak  signal  is 
shown  in  Figure  S.  Two  laser  beams  of  frequencies  cO]  and  C02  =  (0i  -C2  propagate  through  the 
ntmlinear  medium  and  create  a  moving  density  grating  response  that  preferentially  deflects 
energy  from  the  high-frequency  to  the  Icw-frequency  beam.  The  solid  lines  depict  the  crests  of 
the  field-induced  grating  of  wave  vector  q  K]  -  K2.  Dashed  lines  represent  a  spontaneous 
moving  grating  arising  from  thermal  fluctuations  in  the  medium.  Some  of  these  thermal  gratings 
will  be  in  a  prt^  configuration  to  deflect  energy  from  the  high-frequency  to  the  low-frequency 
wave,  in  a  marrner  that  is  indistinguishable  from  the  signal,  to  give  a  noise  component 
represented  by  the  dashed  arrow  in  the  figure.  For  two  incident  plane  waves,  the  total  field  is  of 
the  form 

E(r,t)  =  ejEiCr)  cos[  Kj  •  r  -  ©jt  +  OjCr)  ] 

+  e2E2(r)cos[K2*r-a)2t  +  e2(r)],  (31) 

where  ej,  ^(r),  and  6j(r)  are  the  unit  polarization  vector,  slowly  varying  amplitude,  and  slowly 
varying  phase  of  wave  j. 

The  signal  and  noise  gratings  depicted  schematically  in  Figure  5  are  quantitatively  described  as 
follows:  For  the  signal  grating  (Ref.  8) 

Ae(r,t)  =  aicos(q  •  r-flt)  +  bisin(q  •  r-Qt)  (32) 


where 


ai  =  (ei*e2)e2EiE2/[l  +  (QT)2] 


(33) 


and 
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Figure  5.  Two-wave  mixing  geometry. 


bi  =  -(ei  •  €2)62  EjEj  QtAI  +  (Ox^  (34) 

are  tile  io-phase  and  n/2  out-of-phase  grating  coefficients,  respectively,  for  a  Debye  relaxation 

i 

medium..  For  the  noise  grating,  the  coefficients  a^  and  b(2  are  statistically  distributed  ^cmxling 
to  the  bottom  set  of  equations  in  Table  2  (Ref.  8). 

The  preceding  expressions  may  be  combined  to  obtain  the  total  variation  in  dielectric  constant 
from  both  signal  and  noise  contributions.  When  the  result  is  substituted  into  the  wave  equation, 
the  various  terms  in  the  displacement  vector  give  rise  to  scattering  of  electromagnetic  energy  into 
different  modes,  characterized  by  specific  wave  vectors  and  frequencies.  The  most  significant 
ctmtribution  results  from  those  terms  that  have  minimum  phase  mismatch.  Denoting  these  by 
D(r,t),  (»ie  rinds  (Ref.  8) 

D(r,t)  =  eo  E(r,t)  +  (1/2)  (aj  +  a^)  Ei(r)  cos  ( K2  •  r  -  ©21  +  62)  ej 
+  (1/2)  (ai  +  a^)  E2(r)  cos  ( K 1  •  r  -  ©it  +  0i)  62 
-  (1/2)  (bi  +  bfl)  Ei(r)  sin  ( K2  •  r  -  ©2t  +  62)  C] 

+  (1/2)  (bi  +  bjj)  E2(r)  sin  ( Ki  •  r  -  ©it  +  0i)  62-  (35) 
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Table  2.  Signal  and  noise  gratings  for  two-wave  mixing. 


q  s  Ki  -  K2  fl  =  0)1  -  0)2 

Signal  grating  (field  induced) 

d£(q,Q)  =  aiCOs(q  •  r-Qt)  +  bisin(q  •  r-flt) 

ai  =  62  (ei  •  62)  E1E2  /  [  1  + 

bi  =  -62  (ci  •  C2)  E1E2  flix/Il  +  (fix)2] 

Noise  grating  (thermal! 

aQCOs(q  •  r  -  Qt)  +  basin(q  •  r  -  Qt) 

<  aQ>  =  <  bfl|>  =  0 

<  a^Q  >  =  <  b^Q  >  =  (t/it)  <  i8e(q)|2  >  /  [  1+  (Qt)^] 

=  (8kTte2A^s)/[l+(«t)2] 


Equation  35  describes  an  instantaneous  polarization,  or  snapshot  of  the  medium,  that  includes 
contributions  from  both  signal  and  phase-matched  noise  gratings.  According  to  the  adiabatic 
iqtproximation,  this  polarization  will  remain  sensibly  constant  over  time  periods  that  are  short  in 
ctMoparison  with  the  medium  response  time,  x  =  F'l,  but  still  long  compared  with  the  transit  time 
of  light  through  the  medium. 

Inserting  Equation  35  into  the  wave  equation  (Eq.  16),  making  the  slowly  varying  amplitude  and 
phase  approximation,  and  equating  the  in-phase  and  out-of-phase  terms,  one  finds 

(Ki  •  V)  El  =  k2/4  go(bi  +  bQ)  (ei  •  62)  E2  -  (ao/2)  Ei  (36a) 

(Ki  •  V)  01  =  k2/4  60  (ai  +  ajj)  (ei  •  62)  E2/E1  (36b) 

for  the  pump  wave  and 

(K2  •  V)  E2  =  -K2/4  60  (bi  +  bft)  (ei  •  62)  El  -  (Oo^)  E2  (36c) 
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(Kj  •  V)  62  =  K2/4  So  (aj  +  afl)  (e,  •  ej)  Ej/Ej  (36d) 

for  the  signal.  These  equations  are  identical  to  those  obtained  previously  except  for  the  presence 
of  the  new  noise  tenns,  which  result  in  fluctuations  in  both  amplitude  and  phase.  In  the  case  that 
the  punq)  is  much  stronger  than  the  signal  (E]  »  E2),  these  fluctuations  are  much  more 
important  for  the  signal  than  for  the  pump  beam.  The  last  terms  in  Equations  36a  and  36c  include 
the  effect  of  nonsaturable  background  loss,  where  Oq  is  the  attenuation  coefficient  for  loss  of 
light  intensity.  For  a  nonabsorbing,  isotropic  medium  in  which  losses  are  due  to  scattering  alone, 

oto  =  (034/6jt  c4)  Vs  <  I8e(q)|2  >  =  4/3  (2ic  v/c)^  kT  ej  (37) 

where  v  is  the  optical  frequency. 

Finally,  one  may  obtain  an  expression  for  the  power  dissipated  in  the  medium  during  two-wave 
mixing.  For  a  uniformly  moving  grating,  the  stored  energy  in  the  medium  is  constant  and  the 
dissipated  power  per  unit  volume  is  (Ref.  8) 

P  =  <  E  •  dP/dt>  =  -(1/16JC)  bi(ei  •  e2)EiE2  Q.  (38) 

Equation  38  predicts  an  asymptotic  level  of  power  dissipation  for  Ox  approaching  infinity.  This 
follows  from  the  functional  form  of  b^,  exhibited  in  Equation  34,  which  in  turn  is  a  consequence 
of  the  D4bye  relaxation  model.  As  expected,  the  fluctuating  noise  gratings  make  no  contribution 
to  power  dissipation. 

5.2  EXPERIMENIAL  BESULTS 

Measurements  of  gain  and  noise  were  made  using  liquid  suspensions  of  polytetrafluoroethylene 
(FIFE)  shaped  microparticles  as  artificial  Kerr  media.  Light-induced,  particle-orientation 
rearrangements  account  for  optically-induced  birefringence  with  consequent  large  nonlinear 
optical  response.  Several  studies  have  been  carried  out  to  characterize  the  linear  and  nonlinear 
optical  properties  of  these  suspensions,  and  good  agreement  between  theoretical  predictions  and 
experimental  results  was  found.  This  agreement  is  due  in  large  part  to  the  characteristics  of  the 
suspensions,  which  fit  an  independent  particle-single  orientational  relaxation  time  Debye  model 
remarkably  well.  Particles  are  ellipsoidal  in  shape  (the  dimensions  are  0.4  x  0.2  x  0.2  |xm)  and 
the  particle  polarizability  tensor  has  corresponding  symmetry.  The  particle  suspensions  are 
highly  monodisperse,  and  interparticle  interactions  are  negligible  even  at  relatively  high  volume 
fractions  (0.1-2%).  In  the  present  experiments,  a  1  percent  volume  fraction  of  FIFE  particles 
suspended  in  a  65/35  water/glycerol  host-liquid  mixture  was  typically  used.  This  mixture  was 
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chosen  because  it  matches  the  average  value  of  the  particle  refractive  index  uq  =  1.376,  thus 
minimizing  the  importance  of  light-scattering  fluctuations  due  to  the  isotropic  part  of  the  particle 
polarizability. 

Figure  6  shows  the  experimental  setup  for  nondegenerate  two-wave  mixing.  The  pump  and  the 
probe  waves  are  provided  by  splitting  the  514-nm  TEMqo  output  of  a  cw  Argon-ion  laser  into  a 
1-W  power  pump  beam  and  a  50-mW  probe  beam.  The  two  beams  eventually  intersect  at  the 
focus  of  the  lens  (L)  within  the  sample  at  an  angle  of  5  deg.  The  beam  waist  size  is  about  100  |xm 
and  the  interaction  length  is  0.1  cm.  Before  entering  the  sample,  the  polarization  direction  of  the 
pump  beam  is  rotated  perpendicular  to  that  of  the  probe  so  that  only  the  orientational  particle 
gratings  are  generated  in  the  suspension  by  the  beam  interference  pattern. 


Time  Domain  and 

- ► 

Frequency  Domain 

Data  Analysis 

Laser 

Ar-ion 
514  nm 
cw 


Hgure  6.  Schematic  diagram  of  the  two-wave  mixing  optical  arrangement  and  electronics  for 
noise  analysis.  (L:  focusing  lens;  Ml:  100%  mirror;  PZT:  piezoelectric  transducer, 
D:  photodiode.) 
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The  angular  frequency  shift  Q  required  for  the  nondegenerate  beam  configuration  is.  provided  to 
the  pun^  beam  by  the  mirror  Ml  that  is  mounted  to  a  piezoelectric  transducer  (PZT).  A 
triangular  8-Hz  high-voltage  wave  is  fed  through  the  P2T  to  produce  a  periodic  linear 
displacement  of  the  mirror  with  constant,  sign-inverting  speed.  The  consequent  square-wave-like 
angular  frequency  shift  can  be  varied  by  either  changing  the  amplitude  or  the  frequency  of  the 
mirror  displacement.  The  periodic  intensity  gain  of  the  probe  beam  is  revealed  by  a  low-noise 
photodiode  OD).  A  polarizer  is  placed  before  the  photodiode  to  remove  the  polarized  component 
of  the  pump  light  scattered  into  the  probe  beam  by  the  isotropic  particle  fluctuations.  After  a  low- 
pass  electrical  amplification  (-3  dB  at  3  kHz),  the  photodiode  output  is  fed  both  into  a  digital 
oscilloscope  and  into  a  12-bit,  PC-interfaced  spectrum  analyzer.  For  each  value  of  the  frequency 
shift,  the  average  value  of  the  intensity  gain  was  obtained  by  a  digital  average  on  the  oscilloscope 
while  the  spectral  properties  of  the  gain  fluctuations  were  monitored  on  the  spectrum  analyzer. 
Implementing  this  setup,  particular  care  had  to  be  taken  to  damp  mechanical  vibrations  that  could 
add  random  phase  shifts  to  the  beams  thus  introducing  spurious  noise  in  the  gain  process. 

Figure  7  (top)  shows  the  experimental  results  for  the  probe  intensity  gain  and  noise  fluctuations. 
Pump  light  intensity  in  the  sample  was  approximately  10  kW/cm^  while  the  transmitted  probe 
power  in  the  absence  of  the  pump  beam  was  50  mW.  The  vertical  bars  report  the  transmitted 
probe-beam  power  in  watts  as  a  function  of  the  nondimensional  frequency  shift  Qx,  where 
X  =  Xk  =  18.9  ms  is  the  orientational  response  time  as  determined  through  independent 
measurements.  Each  vertical  bar  is  centered  on  the  average  value  of  the  transmitted  power  and 
the  bar  length  reports  the  root-mean-square  (RMS)  value  of  the  noise  fluctuations  on  the  same 
scale.  The  RMS  noise  fluctuations  are  also  shown  in  the  figure  on  an  expanded  scale  for  clarity 
(full  circles)  and  were  obtained  by  integration  of  the  noise  power  spectrum  over  a  collection 
bandwidth  of  60  Hz.  The  smooth  curve  is  a  fit  to  the  data  using  a(P2)  =  C  [1  +  (Qx)^]'*^,  which 
is  the  form  expected  from  the  theory  and  will  now  be  described. 

For  two  incident  plane  waves,  the  total  field  is  of  the  form  given  by  Equation  31.  Insertion  into 
the  wave  equation,  making  the  slowly  varying  amplitude  and  phase  approximation,  and  equating 
the  in-phase  and  out-of-phase  terms,  gives  a  set  of  four  coupled  equations  (Eq.  36)  describing  the 
evolution  of  the  two  amplitudes  and  two  phases,  which  may  be  integrated  using  the  stochastic 
noise  model  (Ref.  6).  In  the  absence  of  pump  depletion  and  background  loss,  an  analytic  solution 
to  these  equations  may  be  obtained,  which  is  sufficient  for  comparison  with  the  present 
measurements.  Here,  the  solution  is  required  for  the  weak  probe  beam  E2(L),  where  L  is  the 
beam  interaction  length,  in  the  presence  of  a  nondepleted  pump  E].  The  result,  from  Reference  8, 
is 
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Figure  7.  Comparison  of  experimental  and  simulation  results  for  two-wave  mixing  gain  and 
light-scattering  noise. 


E2(L)  =  exp(asL)E2(0)  +  (0^/05  )[exp(asL)  - 1]  Ej 


(39) 

(40) 


and 


aN*(K/4eo)bQ  (41) 

In  Equation  40,  K  =  I  Kjl  =  IK2 1,  Cq  is  the  background  dielectric  constant,  is  the  Kerr 
coefficient,  and  Q  =  O]  -  0)2.  The  last  equality  in  Equation  41  defines  h.  In  Equation  41,  is  the 
quadrature  component  for  a  thermal  fluctuation  grating  having  the  same  orientational 
ccmfiguration  as  the  signal  grating,  thereby  giving  rise  to  a  scattered  light  component  that  is 
indistinguishable  from  the  signal  beam.  Note  that  £2  and  b£2  are  off-diagonal  components  of 
second-rank  tensors,  but  may  be  treated  here  as  scalars  since  in  each  case  only  one  tensor 
component  is  selected  by  the  fixed  polarizations  of  the  writing  beams. 

Inspection  of  Equations  39-41  reveals  that  fluctuations  in  b^  cause  fluctuations  in  the  amplitude 
of  the  transmitted  probe.  In  addition  to  b^,  there  is  an  in-phase  noise  component  a^,  which 
results  in  fluctuations  in  the  phase.  The  mean  and  variances  for  the  grating  amplitude  fluctuations 


(Table  2)  are  derived  in  References: 

<  aQ>  =  <  bQ>  =  0 

(42a) 

<  >  =  (2xB)  <  I5£(q)|2  >  /  [  1  -i-  (Qx)2] 

(42b) 

where  the  dependence  on  frequency  difference  Q  mid  response  time  x  follows  our  assumption  of 
a  single  relaxation  time  Debye  medium.  The  right-hand  side  of  Equation  42b  gives  the  noise 
power  over  a  frequency  bandwidth  B  and  is  obtained  by  multiplying  the  result  of  Reference  8, 
which  gives  the  noise  power  on  a  per  unit  angular  frequency  basis,  by  2kB.  The  right-hand  side 
of  Equation  42b  may  be  evaluated  in  terms  of  the  nonlinear  dielectric  constant,  or  Kerr 
coefficient,  £2  using  Equation  30,  where  q  =  K{  -  K2  is  the  wave  vector  of  the  matched  thomal 
fluctuation  grating  and  Vs  is  the  scattering  or  beam  interaction  volume  determined  from  the 
product  of  the  beam  cross  section  area  A  and  the  interaction  length  L. 

Fot  comparison  with  experiment.  Equation  39  was  evaluated  numerically  using  the  stochastic 
noise  model  (Ref.  6).  The  one  adjustable  paratiKter  in  the  theory  is  the  Kerr  coefficient,  and  this 
is  obtained  here  from  a  fit  to  the  gain  curve,  not  to  the  noise.  The  gain  curve  shown  in  Figure  7 
(bottcnn)  was  obtained  using  £2  =  3.6  x  10~^  cmVerg,  in  reasonable  agreement  with  the  estimate 
for  this  quantity  obtained  from  previous  nonlinear  optics  measurements:  £2  =  2.5  x  lO**^  cmVerg. 
The  remaining  conditions  used  in  the  calculations  are  the  same  as  those  described  above  for  the 
experiment  The  simulation  proceeds  by  the  following  steps:  (1)  A  value  for  Dx  is  selected  and 
Equation  42  is  used  to  obtain  the  mean  and  variance  for  the  noise  amplitude  fluctuations  over  the 
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measurement  bandwidth  B  s  60  Hz.  The  fluctuations  are  assumed  to  be  Gaussian  and  sampling 
is  achieved  using  a  standard  computer  algorithm  incorporating  the  Box-Muller  transformation  for 
die  generation  of  normal  deviates  from  random  numbers  sampled  uniformly  on  die  interval  (0,1). 
(2)  Equation  39  is  evaluated  for  each  sampled  b^  to  give  a  solution  for  E2(L)  and  corresponding 
amplified  signal  power  P2(L).  Steps  1  and  2  are  repeated  on  the  order  of  a  thousand  times  to 
oUain  a  good  statistical  sampling  of  the  noise  over  the  full  Qx  range  of  interest.  Results  of  the 
calculation  are  shown  in  Figure  7  (bottom).  Note  that  both  the  experimental  and  theoredcal  gain 
profiles  follow  frequency  dispersion  curves  characterisdc  of  the  Debye  relaxadon  model.  More 
importantly,  thoe  is  excellent  agreement  between  theory  and  experiment  in  regard  to  the  RMS 
values  of  the  noise  fluctuations,  including  their  dependence  on  fiequency  shift.  This  agreement  is 
remarkable,  in  view  of  the  lack  of  adjustable  parameters  in  the  theory. 

Squaring  the  right-hand  side  of  Equation  39  gives  noise  terms  proportional  to  and  to  In 
most  cases  of  interest,  E2(0)  is  sufficiently  large  that  the  latter  term  may  be  neglected,  leaving  the 
cross  term  containing  the  product  EiE2(0).  Then  one  finds  for  the  standard  deviation  of  the 
power  fluctuations  in  the  amplified  signal  in  the  limit  of  small  gain: 

a(P2)  =  (4jrtB)l/2(hL)l/2  [kT  v  p2(0)]l/2  [1  +  (nx)2]-l/2  (hL  «  1)  (43) 

Since  the  maximum  power  gain  (G  =  e**!-)  in  the  experiment  is  1.05,  the  small  gain 
tqiproximation  used  to  derive  Equation  7  is  valid  for  the  present  discussion.  Equation  43  shows 
the  origin  of  the  [1  +  (Qx)2]-l/2  dependence  seen  in  the  measurements  and  provides  an  explicit 
form  for  the  prefactor  C.  Here  P2(0)  is  the  incident  signal  power,  P2  is  the  fluctuating  signal 
power  at  z  =  L,  and  o(P2)  is  the  standard  deviation  of  P2:  a(P2)  =  1<(P2)^>  -  <P2>^]*^. 

The  solid  curve  in  Figure  7  (bottom)  displays  the  analytical  result  given  by  Equation  43.  The 
open  circles  are  numerical  results  obtained  by  Gaussian  sampling  using  the  stochastic  noise 
model — each  circle  represents  an  RMS  value  computed  for  100  samples  at  a  fixed  value  of  £2x. 
Hie  simulated  noise  is  in  excellent  agreement  with  both  the  analytical  and  experimental  results. 
Finally,  note  that  the  noise  fluctuations  decrease  monotonically  with  frequency  unlike  the 
average  gain,  which  peaks  at  =  1. 

Equation  43  implies  that  the  RMS  amplitude  of  the  noise  fluctuations  increases  with  the  square 
toot  of  die  incittent  probe  power  P2(0).  Figure  8  reports  the  experimental  values  of  the  RMS 
noise  power  (triangles)  as  a  function  of  the  probe  power  for  =  0.  As  can  be  seen,  the  square- 
root  dqiendence  (full  line)  is  well  verified  over  a  decade  range. 
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Figure  8.  The  RN^  noise  power  as  a  function  of  probe  power  for  Qt  =  0.  (The  experimental 
data  points  [triangles]  follow  a  square-root  dependence  on  the  probe  power  [solid 
curve].  The  single  data  point  on  the  (mdinate  gives  the  RMS  noise  power  measured  in 
the  absence  of  the  probe  beam.) 

Theory  predicts  that  the  RMS  probe-power  fluctuations  should  satisfy  Equation  43  as  long  as 
light-induced  gratings  are  dominant  with  respect  to  thermal  gratings.  Equation  39  shows  that  for 
an  incident  probe  power  of  zero,  output  fluctuations  will  occur  proportional  to  Note  that  the 
noise  fluctuations  proportional  to  unlike  those  proportional  to  a{^,  do  not  vanish  on  time 
averaging.  Analysis  similar  to  that  used  to  derive  Equation  43  gives 

<  P2(L)  >  =  (tcxB)  (hL)  (kTv)  [  1  +  (Qx)2]-1  (44) 

where  <  P2(L)  >  is  the  average  power  of  the  noise  at  z  =  L.  The  factor  (kTv)  is  the  Nyquist 
expression  for  the  thermal  noise  power  radiated  by  a  channel  having  a  bandwidth  equal  to  the 
optical  frequency  v.  Equation  44  yields  a  predicted  average  noise  power  at  zero  probe  power  of 
0.2  [iW  due  entirely  to  scattered  pump  light.  The  stochastic  noise  model  sampling  confirms  this 
result  and  gives  an  RMS  power  of  0.28  pW.  The  RMS  noise  power  was  measured  at  zero  signal 
power  by  blocking  the  probe.  The  result  of  this  measurement  (1.8  pW  RMS)  is  represented  by 
the  single  data  point  appearing  on  the  ordinate  of  Figure  8.  At  this  low  noise  level,  the 
discrepancy  between  theory  and  experiment  is  much  greater  than  for  the  previous  comparisons 
made  with  the  probe  beam  on  and  most  likely  is  due  to  extraneous  sources  of  noise  being  present 
in  the  measurement 
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In  summary,  the  first  measurements  of  light-scattering  noise  during  two- wave  mixing  in  a  Kerr 
medium  were  described.  Excellent  quantitative  agreement  between  the  experimental  results  and 
theoretical  predictions  was  obtained,  with  no  adjustable  parameters,  using  the  stochastic  noise 
model  for  thermal  light-scattering  noise  (Ref.  6)  and  its  extension  to  two-wave  mixing  (Ref.  8). 
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6.0  RESEARCH  DIRECTIONS 


Table  3  and  Figures  9-11  describe  the  recendy  completed  results  obtained  for  resonant  media. 
The  model  used  is  that  of  the  damped  harmonic  oscillator.  Noise  arises  from  both  thermal  and 
quantum  fluctuations  of  the  oscillator  polarizability  that  give  rise  to  scattered  light.  Our  previous 
results  for  general  Kerr  media  are  recovered  as  a  special  case  in  the  high  damping  limit.  For 
example,  the  response  curves  sho>vn  in  Figure  1 1  are  identical  to  Equations  33  and  34  for  the 
Debye  relaxation  model  Kerr  medium. 


Table  3.  Light-scattering  noise  in  resonant  media. 

•  Requires  a  unified  description  of  quantum-thermal  noise 

Quantum  coircction 

Se(q.Q)  =  Vs  J  <8e(q.0)  5e(q,t)  >  dt 
Se(q.-Q)  =  Vs  J  <6e(q,t)  8e(q,0)  >  dt 

since  <5p(q,t)  Se(q,0)  >  =  <5e(q,0)  5e(q,t)  >  *  . 

Condition  of  detailed  balance; 

Se(q,-Q)  =  exp  (-ftQ/kT)  Se(q,Q) 

Qassically,  the  correlation  function  is  real  and 
Se(q,-Q)  =  Sc(q,Q). 


Table  3  shows  the  fundamental  nature  of  the  changes  that  need  to  be  made  in  the  fluctuation 
theorem  to  include  quantum  effects.  The  results  described  thus  far  in  this  repon  apply  to  the 
thermal  noise  limit  (Q  «  kT/A)  typical  of  slow  media  having  large  nonlinear  optical  coefficients. 
In  this  linoit,  the  maximum  fi^uency  difference  Q  is  limited  by  the  medium  response  time  and 
quantum  ncnse  effects,  which  can  be  dominant  for  Raman  amplifiers  (Ref.  1 1),  may  be  neglected 


29 


Results  for  Beta  »  0.1;  wO  •  1 
Top:  Gain  Curve 


Figure  9.  Results  for  damped  oscillator,  low  damping. 

as  these  are  oiders  of  magnitude  smaller  than  the  noise  described  here.  A  soon-to-be-completed 
imblicadon  will  extend  dw  theory  to  cases  for  which  Q  is  comparable  to  or  greater  than  kT/h  in  a 
unified  treatment  of  quantum-thermal  noise. 

Widi  further  development,  the  results  may  serve  as  the  scientific  basis  necessary  for  the  design  of 
new  low-noise  nonlinear  optical  materials — ^probably  with  best  prospects  in  the  general  classes  of 
phottxefractive  polymers  and  charge  transfer  media.  Furthermore,  the  results  can  provide  the 
foundation  for  investigations  into  the  noise  properties  of  optical  communication  channels  and 
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Figure  10.  Results  for  damped  oscillator,  intermediate  damping. 

optical  storage  media.  These  investigations  can  include  determination  of  the  Shannon  channel 
capacity  (some  woric  along  these  lines  has  already  been  carried  out  for  photorefractive  media) 
and  new  aj^lications  of  c^dcal  signal  processing  and  computing  (e.g.,  the  optical  wavelet 
transform  and  t^tical  methods  for  signal  and  image  compression). 


Results  for  Beta  •  3.1; 


Hgure  1 1.  Results  for  damped  oscillator,  high  damping. 
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